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04 : Abstract 
C . 

! & . We derive a general formula for the center-of-mass (CM) energy for the near-horizon collision 

\ of two particles of the same rest mass on the equatorial plane around a Kerr black hole. We then 

q ' apply this formula to a particle which plunges from the innermost stable circular orbit (ISCO) and 

collides with another particle near the horizon. It is found that the maximum value of the CM 
energy E cm is given by E cm / (2mo) ~ 1.40/ y 1 — a 2 for a nearly maximally rotating black hole, 
where mo is the rest mass of each particle and a* is the nondimensional Kerr parameter. This 
■ coincides with the known upper bound for a particle which begins at rest at infinity within a factor 

. of 2. Moreover, we also consider the collision of a particle orbiting the ISCO with another particle 



X 



on the ISCO and find that the maximum CM energy is then given by E cm / (2mo) — 1.77/ \J 1 — a\ . 
In view of the astrophysical significance of the ISCO, this result implies that particles can collide 
around a rotating black hole with an arbitrarily high CM energy without any artificial fine-tuning 
in an astrophysical context if we can take the maximal limit of the black hole spin or a* — > 1. 
On the other hand, even if we take Thome's bound on the spin parameter into account, highly or 
moderately relativistic collisions are expected to occur quite naturally, for E cm / (2mo) takes 6.95 
(maximum) and 3.86 (generic) near the horizon and 4.11 (maximum) and 2.43 (generic) on the 
ISCO for a* = 0.998. This implies that high-velocity collisions of compact objects are naturally 
expected around a rapidly rotating supermassive black hole. Implications to accretion flows onto 
a rapidly rotating black hole are also discussed. 
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I. INTRODUCTION 



Recently, Banados, Silk and West [lj] showed that if two particles which begin at rest 
at infinity collide near the horizon of a maximally rotating Kerr black hole ^ and if the 
angular momentum of either particle is fine-tuned, the center-of-mass (CM) energy E cm 
of the two particles can be arbitrarily high and hence the maximally rotating black hole 
might be regarded as a Planck-energy-scale collider. We here call this the Banados-Silk- 
West (BSW) effect. This scenario was subsequently criticized from several points, 

such as astrophysical bounds on the black hole spin parameter, the effects of gravitational 
waves, the self-gravity of the particles and the long proper time needed for such a collision. 
In the near-maximal rotation, the CM energy of two particles of mass m is bounded by 
-E'cm/(2mo) ~ 2.41/ ^/l — a%, where a* is the nondimensional Kerr parameter 3|, |4|. On the 
other hand, Grib and Pavlov j^j proposed a slightly different mechanism that E cm can be 
arbitrarily high even in the non-maximal rotation if the particle has experienced multiple 
scattering and fine-tuned its angular momentum before the relevant collision. The BSW 
effect is also analyzed in the Kerr-Newman family of black holes js], general stationary and 
axisymmetric black holes 3|, and static charged black holes 

In the BSW effect, to obtain such an arbitrarily high E cm , the angular momentum of 
either particle must be fine-tuned. However, there is a natural mechanism in astrophysics 
to tune the particle's energy and angular momentum. This is the innermost stable circular 
orbit (ISCO). The ISCO around a Kerr black hole is studied in detail by Bardeen, Press and 
Teukolsky [9]. 



In the geometrically thin and optically thick accretion disk model 10, lH|, which is known 
as the standard accretion disk model, a rotating fluid or plasma gradually takes a circular 
orbit which is closer to the black hole as the fluid transfers its angular momentum outwards 
and releases its energy by electromagnetic radiation in the time scale of viscosity, which 
is much longer than the black hole dynamical time scale. This electromagnetic emission 
can be observed by radio interferometers and X-ray observation satellites. Once the fluid 
reaches the inner edge of the accretion disk, which is given by the ISCO, it begins to plunge 



into the black hole in the dynamical time scale [26]. In the plunging phase, its energy 
and angular momentum are approximately conserved. It should be noted that due to the 
accretion of radiation emitted from the disk, there is an astrophysical upper bound on the 
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nondimensional Kerr parameter, what we call Thome's bound, |a*| < 0.998 [l^ . 

Another important example where the ISCO plays a crucial role is inspirals of stellar- 
mass compact objects into supermassive black holes, which are called extreme mass-ratio 



inspirals. Extreme mass-ratio inspira 
Laser Interferometric Space Antenna 



s are interesting sources of gravitational waves for the 



151 ] . If instead the central mass is an intermediate-mass 



black hole, these are interesting sources for the DECi-hertz 



nterferometer Gravitational 



wave Observatory [16|, [17| and the Big Bang Observatory [18|, [19]. Also in this case, an 
inspiralling compact object gradually takes a circular orbit which is closer to the black 
hole as the object transfers its angular momentum outwards and releases its energy by 
gravitational waves in the time scale of gravitational radiation, which is much longer than 
the dynamical time scale. Once the compact object reaches the ISCO, it begins to plunge 
into the black hole in the dynamical time scale. In the plunging phase, its energy and 
angular momentum are approximately conserved. 

In this paper, we derive a general formula for the CM energy for the near-horizon collision 
of two particles on the equatorial plane around a Kerr black hole, which is valid in both 
the maximal and non-maximal rotation cases. Then, we apply this formula to the near- 
horizon collision of two particles, either of which is plunging from the ISCO. We find that 
the BSW effect occurs in the near-maximal rotation and that the maximum value for the 
CM energy of the ISCO particles is the same as the upper bound for the particles initially at 
rest at infinity within a factor of 2. This implies that the BSW effect is not an artificial but 
physically realistic astrophysical phenomenon. We also consider the collision of a particle 
orbiting the ISCO with another generic particle on the ISCO and find that the associated 
CM energy can also be arbitrarily high in exactly the same sense as BSWs, although the 
dependence on the black hole spin parameter is quite different. We neglect the effects of 
gravitational waves and the self-gravity of the particles. 

This paper is organized as follows. In Sec. II we briefly review particle orbits and the 
CM energy for the collision of two particles in the Kerr spacetime. In Sec. Ill, we discuss 
particle orbits near the horizon and derive a general formula for the CM energy for the 
near-horizon particle collision. In Sec. IV, we apply this formula to a particle which plunges 
from the ISCO and obtain the CM energy for different collisions. In Sec. V, we investigate 
the collision of a particle orbiting the ISCO with another particle on the ISCO. Section VI is 
devoted to conclusion and discussion. We use the units in which c = G = 1 and the abstract 
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index notation of Wald 



20] 



II. CM ENERGY FOR PARTICLE COLLISION IN THE KERR SPACETIME 



In this section, we briefly review particle orbits and the CM energy for the two-particle 
collision on the equatorial plane of the Kerr spacetime in the general situation, following 
4j, |5[. We use a similar notation to that of Grib and Pavlov [5]. 
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A. Particle orbits in the Kerr spacetime 

The line element in the Kerr spacetime in the Boyer-Lindquist coordinates is given by 



Ei] 



2Mr\ , 2 4Mar sin 2 9 , p 2 , 2 2 in2 
ds z = - ( 1 — \ dt 2 — 2 d<j)dt + ^-dr 2 + p 2 d6 2 



P' 



2 2 2Mra 2 sin 2 #\ . 2 2 
( r 2 + a 2 H — 2 J sin 2 6>d0 2 , 

where a and M are, respectively, the spin and mass parameters, p 2 = r 2 + a 2 cos 2 # and 
A = r 2 - 2Mr + a 2 . We assume a > without loss of generality. If a 2 < M 2 , A vanishes 
at r = r± = M ± V^f 2 — a 2 , where r = r + and r = r_ correspond to an event horizon and 
a Cauchy horizon, respectively. Here, we denote r + = r# and r_ = rc- In this coordinate 
system, the time translational and the axial Killing vectors are respectively given by 



A, 

The surface gravity of the Kerr black hole is given by 



K = ^ — 



Thus, the black hole has a vanishing surface gravity and hence is extremal for the maximal 
rotation a 2 = M 2 , while it is subextremal for the non- maximal rotation a 2 < M 2 . The 
angular velocity of the horizon is given by 

a„ a 



r 2 H + a 2 2M(M + \J M 2 — a 2 ) ' 
The Killing vector \ a = £, a + ^h^ 0- is a null generator of the event horizon. 
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Here we concentrate on particles which move on the equatorial plane 9 = tt/2. Then, the 
four velocity u a = x a of the particle has a vanishing 9 component, i.e., u e = 0, where the 
dot denotes the differentiation with respect to the affine parameter of the geodesic. From 
Eq. (]2.ip . the line element on the equatorial plane in the Kerr spacetime is given by 



ds 2 



1 - 



2M 



dt 2 



4aM 



— dtd(f) + —dr 2 + [r 2 + a 2 



2Ma s 



(2.1) 



v I r A \ r 

Associated with the Killing vectors £ a and ip a , we have the following conserved quantities 
along a geodesic on the equatorial plane: 



-9abi a u b = -ut = -{gttu 1 + gnu*), 



g a b^ a u b = — g^u 1 + g^u* 



(2.2) 
(2.3) 



where e and L correspond the specific energy and angular momentum, respectively. Solving 
the above for it* and ir\ we have 



t 



2Ma 2 \ 2Ma 
r 2 + a 2 H ) e L 

r 



r 

2M\ 2Ma 
1 I L + e 



To have t > 0, the condition 



1 

A 
1 

A 



9 9 2Ma 2 \ 2Ma T 
r 2 + a 2 H | e L > 



(2.4) 
(2.5) 



must be satisfied outside the event horizon. In the limit to the horizon r — > from outside, 
this condition reduces to 



/ < / 



2 1 + ,/1-a 



ii 



or 



e > QhL, 



where we put a* = a/M and / = L/M. In terms of a*, < a* < 1 for the subextremal case, 
while a* = 1 for the extremal case. 

Substituting Eqs. f!2.4[) and (12. 5p into the normalization condition u a u a = — 1 with 9 = 
7r/2 and u e = 0, or 



- 1 



2M\ , t 2 AaM t A r 2 , „ n { n n 2Ma 2 



u 



■11V + -M 2 + [r 2 + a 2 + 



-1, 



we obtain 



\r 2 + V cS (r) = 0, (2.6) 



where the effective potential V e g(r) is given by 



M L 2 -a 2 (e 2 -l) M(L - aef e 2 - 1 
Mr) = ~ + ^ 7i —■ (2-7) 



The effective potential can be efficiently analyzed by introducing y = M/r. We put 

g(y) = -2(1 - a*e)V + [I 2 - a\(e 2 - l)]y 2 - 2y - (e 2 - 1) (2.8) 

and 

D{y) = aly 2 -2y+l. 

Then, we have V e s = g(y)/2 and A = r 2 D(y). There are two positive roots of D = for 
< a* < 1, which are given by 

f- y/T^ 2 1 + y/T 

vh = 2 — ' yc = T 



al 



where yu = M/r^ and yc = M/rc correspond to the event horizon and the Cauchy horizon, 
respectively. These two roots coincide with each other at y — 1 in the extremal case a* = f . 
The region outside the horizon is transformed into < y < yu- We should note the following 
useful relation 

For a particle which is initially at rest at infinity, i.e., marginally bound e = 1, to reach 
the horizon, the potential g(y) must be nonpositive for < y < yu- For e = 1, the potential 
is given by 

g{y) = -y[2(l-a,fy 2 -l 2 y + c A- 

Thus, the condition reduces to that 2(1 — a*) 2 y 2 — l 2 y + 2 is nonnegative. After some 
straightforward calculation, we can obtain the following condition 5[] 



-2(1 + VT+aZ) = I l <1<Ir = 2(1 + VT^O- 

A similar condition also exists for the non-marginally bound case. However, we should note 
that this does not apply if the particle scatters with other particles and changes its energy 
and angular momentum on the way to the horizon. 
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B. CM energy of two particles in the Kerr spacetime 



We consider the collision of two particles 1 and 2 of the same rest mass mo- We assume 
that the two particles are at the same spacetime point. The four momentum of particle % 
(i = 1, 2) is given by 

Pi = m ul 

where uf is the four velocity of particle i. The sum of the two momenta is given by 

P a t=p a i+P2- 

The CM energy E cm of the two particles is then given by 

E 2 cm = -ptPta = 2mg(l - g ab u>\)- (2.9) 

On the background metric (j2TT|) . using Eqs. (TjOJ, Q, §ZM) and (EU) in Eq. (|23)» . the 
CM energy of two particles 1 and 2 in the Kerr spacetime is calculated as 
F 2 

era i t t / t <t> i <fr t\ rr 4> <t> 

^—2 = 1 - g tt u x u 2 - gt^u^l + u\u 2 ) - g rr u x u 2 - g^u^ 
F(y)-G(y) 

= 1 ~ 6162 + D(y) ' (2 ' 10) 
where and k are e and / for particle i, 

F(y) = 2[aly 2 (1 + y) + (1 - y)]e x e 2 - 2a*y\e l l 2 + l x e 2 ) - (1 - 2y)y 2 / 1 / 2 , (2.11) 

G = (±y/=n){±y/=gd, (2.12) 

^i(y) = -2{h - a^) 2 y 3 + [l* - a\{e 2 - l)}y 2 - 2y - (e 2 - 1), (2.13) 

D(y) = aly 2 -2y + l, (2.14) 

and the sign in front of \/—gi in the expression of G in Eq. (I2.12p corresponds to the sign 
of u\. In the following we assume f < for both particles and hence G(y) = a/ (— gi)(— g~2)- 

III. NEAR-HORIZON COLLISION AROUND A KERR BLACK HOLE 

A. Near-horizon collision around a subextremal Kerr black hole 

We will see the near-horizon behavior of particles with the angular momentum I = Ijj, 
which we call critical, and smaller angular momentum I < Ih, which we call subcritical. We 
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find 

g(y H ) = -y 2 H {2e - a*y H l) 2 = -aly H {l H - l)\ (3.1) 

noting that Ih can be written in terms of a*, yn and e as Ih = ^ejia^yn)- Hence, giyn) < 0. 
It is interesting to see whether a particle with I = Ih which approaches the horizon is 
possible. With I = Ih, we have g(yn) = 0, while 



g'{y H ) = -2^— ^K 1 + V 7 ^!) V + «*], (3.2) 
a * 

where the prime denotes the differentiation with respect to the argument. This is negative 
for the subextremal Kerr case a 2 < 1. This means that for the subextremal case, the effective 
potential g(y) is positive in the vicinity of the horizon and hence a particle with the angular 
momentum I = Ih is prohibited to approach the horizon. On the other hand, there does 
exist a particle with slightly smaller angular momentum I = Ih — 5 which approaches the 
horizon in the vicinity of the horizon. The CM energy for the collision involving this particle 
can be arbitrarily high in the limit 5 — > even in the subextremal Kerr case 5[. For the 
subcritical orbit / < Ih, Eq. (13.1 ft implies that r is given near the horizon in terms of the 
particle's proper time r as 

r — r H — —a*y H (l H — l)r + const. (3.3) 

This means that for the fixed initial radius, a subcritical particle reaches the horizon after 
a proper time inversely proportional to (Ih — I)- 

Then, we will take the limit to the horizon in Eq. (I2.10p to consider the collision near the 
horizon. Noting 

F(y H ) = aly H ( l m- hWm-h), 

where lm is the critical angular momentum Ih for particle i ii = 1,2), combined with 
Eq. ( 13.11) . we can see that the terms of 0(1) in the numerator F — G of the fraction on 
the right-hand side of Eq. (12.101) cancel out. The nonvanishing contribution comes from the 
next order terms. Using rHospital's rule, the result is the following: 

EL . F'-G' 



cm 



1 — + lim 



2ml y ^ VH D' 
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The derivatives are given by 



F'(y) = 2[al(2y + 3y 2 ) - l] ei e 2 - 6a^ 2 ( ei Z 2 + l±e 2 ) - 2y{\ - 3y)hl 2 , 

g[{y) = -6(h - a^fy 2 + 2[l 2 - al(e 2 - l)]y - 2, 

D\y) = 2{aly-l), 

G' = 



2 \gi 92, 

From this form, we can see that there are two first-order poles, where g%{yii) = for i = 1,2. 
By implementing the calculation and taking the limit, we reach the following formula: 



E cm 4[(l m - ZQ - (l H 2 - k)] 2 + {}mh - Imh) 2 . , 

2m V + 16(lm-h)(l H 2-h) ' 1 ' ' 

This is the formula for the CM energy of two particles along the general geodesic orbits on 
the equatorial plane. We should note that the right-hand side is given only in terms of the 
particles' angular momenta l\ and l 2 and their critical values I hi and Z# 2 . in terms of the 
quantities which have more direct physical meanings, Eq. ( 13. 4 p can be rewritten as follows: 



E cm 4M 2 [( ei - fjglq) - (eg - £l H L 2 )f + ( ei L 2 - e 2 LQ 2 

2m V 16M 2 ( ei -^i)(e 2 -^L 2 ) ' l ' J 

where Lj is L for particle z. In fact, as we will prove in Sec. IIIIBl Eq. (13. 4p or equivalently 
Eq. (13.51) is valid even for the extremal Kerr black hole simply by taking the near-extremal 
limit a* — > 1. The necessary condition for obtaining an arbitrarily high E cm is therefore 
I — > In or VLhL — > e for either of the two particles. 

If we assume that only particle 1 is near-critical in Eq. ( 13. 4(1 . we obtain 



Ecm /4 + Zf^ — h 



2m y 16 — ^i 
For ei = e 2 = e, we denote Ihi = Ih2 = In and Eq. (13.41) reduces to 



(3.6) 



£ cm _ i (h-i 2 m+i 2 H ) 



(3.7) 



2m V 16(lH-h)(l H -l2) 

which reproduces the corresponding formula in jlj]. When we set e\ — e 2 — 1, h — Ir and 
l 2 = II in Eq. (13. 7p . we obtain 



£cm 1 /(l-a2) + (l + vTT^+vT^;)^ 



2^0 ^r^!V i + v^ 
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a? 



(3.8) 



This reproduces the corresponding formula in j^J. This provides an upper bound for the 
collision of two marginally bound particles. For a* = 0.998, E cia / '(2mo) — 9.49 for this 
collision. In the limit a* — > 1, we have 



1 + ^ ^ 2 - 41 f39l 

2m yr^! </r^T 

which are given in |3|-|5|. 

If t\ = 1 and /i = Ir for particle 1 and particle 2 takes a subcritical orbit, we obtain 

^ » i V2 fiZi 2 (3.10) 

in the near-extremal limit. Although the numerical factor depends on the choice of e2 and 
I2, the proportionality to (1 — a 2 )" 1 / 4 does not change as long as e\ = 1 and Zi = Ir. 

B. Near-horizon collision around an extremal Kerr black hole 

For the extremal Kerr black hole a 2 = 1, the effective potential is given by 

g(y) = -2(1 - e)V + (/ 2 - e 2 + l)y 2 - 2y - (e 2 - 1). (3.11) 

The double root y = yn = 1 of D(y) gives an event horizon. We should note that the region 
outside the horizon, r > M, is transformed to < y < 1. 

From Eqs. (13. ip and (13. 2p . for the critical orbit I = In = 2e, we find g(yn) = g'ijju) = 
and hence y = yu is a stationary point of the effective potential. In fact, for the critical 
orbit, we have 

g{y) = -(3e 2 - 1)(1 - y) 2 + 2e 2 (l - y) 3 . (3.12) 

Thus, there exists a critical orbit in the vicinity of the horizon if and only if 3e 2 > 1 and 
then the effective potential takes a maximum which is zero at y — yn — 1- From this fact, 
one might infer an unstable circular orbit for a massive particle at y — yn, i-e., on the 
horizon which is a null hypersurface. This apparent paradox is resolved in Sec. IHI CI On 
the other hand, the maximal point on the horizon implies the existence of an orbit for a 
massive particle with / = In which asymptotes the horizon. In fact, from Eq. (I3.12p we have 



2 e 2 

-v / 3^I(l -y)\\\- 3^37(1 - y). (3.13) 
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and this can be integrated to give 

In \r — M\ ~ — VSe 2 — l-^ + const, 
near the horizon. Thus, the critical particle approaches the horizon as r — > oo, as shown 



in 



5|. For the sub critical orbit / < 2e, since 
-g(y) = [(2e - I) - 2(e - 0(1 - y)f - [1 - (e - /)(3e - Z)] (1 - y) 2 - 2(e - /) 2 (1 - yf 
we have a different behavior of u r as 



r = W = - [(2e - I) - 2(e - J)(l - »)] ,/l - 11 = (e = = + 2(<? = = "'(1 - »)». 

[(2e-Z)-2(e-Z)(l-j/)] 2 



The proper time for the subcritical particle to reach the horizon is inversely proportional 
to (2e — /) because Eq. (13.31) is still valid even in the extremal Kerr case. 

If particle 1 takes a critical orbit but particle 2 takes a subcritical orbit, the CM energy 
is given by the near-horizon limit of Eq. (I2.10p as 



^cm ^ / (2e 2 -/ 2 )(2 ei - y/3^T) 

2m ~ V 2(1 -y) ' 1 ' ] 

where we have used Eq. (13.131) . For the special case e± = e 2 = e, Eq. (I3.14p reproduces the 
corresponding formula in [5]. 

If both particles take subcritical orbits, the fraction on the right-hand side of Eq. f l2.10p 
is bounded because both the numerator F — G and the denominator D have a second-order 
zero at the horizon y = yjj = 1. Estimating the terms of 0((1 — y) 2 ) in F — G and D by 
Taylor series expansion, we obtain 



^cm _ /l 
\ 2 



2e 2 - l 2 1 + e\ 2ei - h 1 + e| 

1 — eie 2 H 1 

2ei - h 2 2e 2 - 1 2 2 



(3.15) 



More systematically, we can take the following approach. For the extremal case, the numer- 
ator F — G and the denominator D both must have a second-order zero at y = yjj. Using 
l'Hospital's rule twice, we obtain 



rp2 Z7>// ri» 

_2| = 1 - ei e 2 + lim — , 3.16 
2mo y^^H -D" 
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where the second-order derivatives are given by 

F"(y) = 4[a 2 (l + 3y)] ei e 2 - Ua^yfak + he 2 ) - 2(1 - 6y)hl 2 , 
g'!{y) = -I2{l l -a,e l fy + 2[/ 2 - a 2 (e 2 - 1)], 
D"{y) = 2al 



G" = G 





(A 




\ 


(A 






U 


92 J 




\9i 


92 J 



and a* = 1 for the extremal case. It is found that this approach also yields Eq. f 13 . 1 5 j) . 

In the course of derivation, it is not so obvious whether the formula for the subextremal 
Kerr black hole given by Eq. (13. 4p reproduces the formula (13. 15[) for the extremal case if we 
take the near-extremal limit a* — > 1 in the former. In fact, it is not difficult to see that this 
is the case by putting lm = 2e; in Eq. (13.4R . Therefore, the general formula (13.41) or (13. 5ft . 
which has been derived for the subextremal case, is applicable in both the subextremal and 
extremal cases. 

We can confirm that for the special case e\ = e 2 = e, Eq. ( 13. 1 5j) reduces to 

E cm j 1 + e 2 (h-l 2 ) 2 
2m y 4 (26-/0(26-/2)' 

which coincides with the corresponding formula in Moreover, Eq. (13. 15[) reduces to 

£cm j I ( 2-h 2-/ 2 \ 

2m \j2\2-k 2-/J' 
for the special case e x = e 2 — 1, which reproduces the formula discovered by BSW [l|. 



C. The circular timelike orbit on the extremal Kerr black hole horizon is fake 

As we have seen in Sec. IIIIBl in the extremal Kerr case, r = r# is a zero and maximal 
point of the effective potential. From this fact, one might infer that an unstable circular 
orbit for a massive particle is possible at r = th- However, we will show that this is not 
real. 

We should note that the Boyer-Lindquist coordinate system has a coordinate singularity 
at r = th- To avoid the complication due to the coordinate singularity, we move to the 



ingoing Kerr coordinates 211 ] : 



dv dv 
dv = dt + (r 2 + a2 )"^"? dip = dcf) + a—. 
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The line element then can be written as 

j2 / 2Mr\ , 2 OJ ^ 2J/J2 [(r 2 + a 2 ) 2 -a 2 Asin 2 #]sin 2 tf 2 
ds 2 = — ( 1 — cfo 2 + 2dudr + p 2 d6» 2 + ^ - s - d(p 2 

V p J p 2 



—2a sin 2 Odpdr — — sin 2 Odcpdv. 



AaMr 
P 

The Killing vectors are given by 



On the equatorial plane 9 = tt/2, the line element in the extremal Kerr spacetime is given 
by 

/ 2M\ ( 2M 3 \ AM 2 
ds 2 = - 1 )dv 2 + 2cfodr + r 2 + M 2 + rfy? 2 - 2Mdipdr dipdv. 



r 

The conserved quantities are given by 



e=-^=(l-^)*-^ + ^, (3.17) 

/ 2M 3 \ 2M 2 
L = g ah i\) a u h = I r 2 + M 2 + — — J - Mr —v. (3.18) 

Putting r = th = M in Eqs. (13. 17)) and ( 13 . 1 Sj) . we have 

e = _ t ;_r + 2Mv9, (3.19) 
L = -2Mu - Mr + AM 2 ip. (3.20) 

The norm of u a can be written at r = M as 

u a u a = v 2 + 2rv + 4M 2 ip 2 - 2M(pr - 4M(pv. (3.21) 

If we assume that the particle remains on the horizon r = M, we have r = and we may 
conclude L/e = 2M or e = L = from Eqs. (13.191) and (13. 20[) . On the other hand, one 
cannot solve Eqs. (13. 19[) and (13.201) for v and ip separately in terms of e and L because of 
the degeneracy and hence cannot obtain Eqs. ( 12. 6 p and (12. 7[) . In other words, Eqs. ( 12 .6p and 
( 12. 7p do not make sense in the present case. Instead, we obtain from Eqs. ( I3.19P and ( I3.2ip 

u a u a = e 2 . 

Since u a u a < for causal geodesies, this means that a causal geodesic can remain on the 
horizon only if it is a null geodesic with e = L = 0. Thus, the timelike circular orbit on 
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the horizon, which might be inferred from the stationary point of the effective potential at 
r = Thi is fake. Note also that the angular velocity of this null geodesic which remains on 
the horizon is given by 

dip = Lp 1 _ n 
dv v 2M H ' 
indicating that this null geodesic is a generator of the event horizon. 

This is of course entirely consistent with the following general argument. Any timelike 
curve cannot remain on an event horizon because the event horizon is normal to the Killing 
vector x a i which is null on the horizon, and hence the tangent space at a point on the horizon 
is spanned by \ a an d two spacelike vectors and e" 2 ) which are orthogonal to x a - Any 
linear combination of the null vector \ a an d the spacelike vectors e%\ and e@) is either null 
or spacelike. 



IV. NEAR-HORIZON COLLISION OF A PARTICLE PLUNGING FROM THE 
ISCO 



As is emphasized in |5| and in Sec. Ill A[ the upper bound (I3.9P applies only for the 
particles which begin at rest at infinity and reach the horizon all the way from infinity along 
the geodesic. If a particle loses or gains its energy or angular momentum on the way to the 
horizon, this limit does not apply. In this respect, particles plunging from the ISCO are 
considered very natural as particles plunging into the horizon in an astrophysical context. 



The ISCO in the Kerr spacetime is explicitly given by Bardeen, Press and Teukolsky {9]. 

The circular orbit on the equatorial plane in the Kerr metric is given by V e g(r) = V^ s (r) = 0, 

where the prime denotes the derivative with respect to the argument. The condition implies 

r x /2( r _ 2M) + saM 1 / 2 
6 ~ r 3 / 4 (r 3 / 2 - ZMr 1 / 2 + s 2aM 1 / 2 ) 1 / 2 ' ^ ' ' 

M x l 2 {r 2 + a 2 — s2M 1 / 2 ar 1 / 2 ) 
~ s r 3/4( r 3/2 _ zMr 1 ' 2 + s2aM 1 /2)i/2' ( ' ' 

where we have assumed < a < M and s — 1 and —1 correspond to the prograde and 
retrograde orbits, respectively. The ISCO is determined by the condition de/dr = dL/dr = 
0. The radius of the ISCO is then given by 



3 + Z 2 - s[(3 - Z 1 )(3 + Z t + 2Z 2 )} 1 ' 2 , (4.3) 



r ISCO 

M 

Z x = 1 + (1 - a^) 1/3 [(l - a*) 1 ^ + (1 + a#) i/ 3]) Z2 = (3a 2 + Z 2 ) i/ 2) (44) 
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where < a* < 1. The energy and angular momentum of the particle at the ISCO are 
calculated by substituting Eqs. (14. 3[) and (14. 4p for r = risco into Eqs. (14. ip and (I4.2p . 

To see the behavior in the near-extremal limit a* — > 1, we put a* = 1 — e and expand the 
above obtained expression in terms of e. From Eqs. (14.41) . we obtain 

Z x = 1 + 2 2 /V/ 3 + 2 1 / 3 e 2 / 3 + 0(e), Z 2 = 2 + -2 2 / 3 e^ + 7 -2^e 2 / 3 + 0(e). 

2 8 

Then, using Eqs. (ETTp . (jOl and ( Ojl . we obtain 

!^CO = j 2 2/3 e l/3 Z 2 l/3 e 2/3 0(f) (45) 

e = ^ + ^VSe 1 / 3 - ^^e 2 / 3 + 0(e), (4.6) 

I = - V3 + V 3 ^ 1 / 3 + i2 1 / 3 v / 3e 2 / 3 + 0(e), (4.7) 
3 3 6 



and hence 



- = 2 + 3-2 1 / 3 e 2 / 3 + 0(e) 
e 



for the prograde orbit (s = 1). For the retrograde orbit (s = —1), we obtain 

^ =9 + 0(6), 

e = ^+0(4 
22 

Z = -yV / 3 + 0(e), 

and hence 

^-f + ow. 

e 5 



On the other hand, Zjj can be written as 

Ih 2a* 



2 + 2v / 2e 1/2 + 0(e). 



e 1 - ^1 - a 2 

Therefore, the prograde ISCO particle has the angular momentum which coincides with the 
critical value Ih in the near-extremal limit. 

We should note that for the marginally bound particle with I = we have l/e = In = 
2 + 2c 1 / 2 . Therefore, (l H - l)/e « 2 v / 2e 1 / 2 for the prograde ISCO, while (l H - l)/e « 
2(\/2 — l)e x / 2 for a marginally bound particle with / = In- Since all other factors are 
nonzero finite, a particle which plunges from the prograde ISCO collides with a generic 
particle with E cm oc e -1 / 4 , as a marginally bound particle with I = l R does. It also follows 
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that if a particle plunging from the prograde ISCO collides with a marginally bound particle 
with / = Ir, the CM energy is bounded even in the near-extremal limit a* — > 1. 

Using Eq. (13. 6p with Eqs. (l4.5p -( l4~Tj) . we can easily estimate the CM energy near the 
horizon for the near-extremal Kerr black hole. If particle 1 is a particle plunging from the 
prograde ISCO and particle 2 takes a subcritical orbit, we obtain 

2m ~ 2V2 3 V4 ■ ^ • > 

Thus, the CM energy can be unboundedly high in the limit a* — > 1. Since the dependence 
(1 — a*) -1 / 4 is common to the upper bound (13. 9p for marginally bound particles, we can 
conclude that the BSW effect occurs for a particle plunging from the prograde ISCO and 
in this case the fine-tuning of the energy and the angular momentum is naturally realized 
in the standard accretion disks with electromagnetic radiation or in inspiralling binaries 
with gravitational wave radiation. The ratio of E cm for the ISCO particle to that for the 
marginally bound particle with I = Ir is given ~ 0.582 in the near-extremal 

limit. 

In the following, we consider the near-horizon collisions of a particle plunging from the 
prograde ISCO with (a) a marginally bound particle with I = II, (b) a particle plunging from 
the retrograde ISCO, and (c) a marginally bound particle with / = 0. In the near-extremal 
limit a* — > 1, we obtain 

2^~ yT^f 

where the numerical factor a = \/2e2 — h/ (2 1 / 2 3 1//4 ) is calculated to be J (2 + y/2) / \/3 ~ 
1.40, 4/3 ~ 1.33 and 1/^3 ~ 0.760 for cases (a), (b), and (c), respectively. 

For the general values of a* in < a* < 1, we can calculate the CM energy using Eq. (13. 4p 
with Eqs. (I4.ip - (l4.4p . The result is summarized in Fig. [TJ where -E cm /(2m ) is multiplied 
by a/1 — clI for clarity. In this figure, the solid, dashed, and dotted curves denote cases (a), 
(b), and (c), respectively. The CM energy for two particles, either of which is a particle 
plunging from the prograde ISCO is always below the upper bound (13. 8p for marginally 
bound particles. We should note that the dependence of (-E cm /(2m ))(l — a*) 1,/4 on a* 
in < a* < 1 is very weak for cases (a), (b), and (c) and hence the empirical formula 
-E-cm/ (2mo) ~ 1/ \J\ — al is a very good approximation within a factor of 2 or so for all 
values of a* in < a* < 1. Thus, this formula provides the typical value for the CM energy 
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for the near-horizon collision of a particle which plunges from the prograde ISCO with a 
generic subcritical particle. For the near-maximal rotation, the maximum value for the CM 
energy coincides with the upper bound for a marginally bound particle within a factor of 2. 

If we use Thome's bound a* = 0.998 for the spin parameter, E cm /(2mo) is calculated to 
be 6.95, 6.61, and 3.86 for cases (a), (b), and (c), respectively. This means that a highly 
relativistic collision can naturally occur near the horizon of a rapidly rotating black hole in 
an astrophysical context. Note that with highly relativistic collision we here mean systems 
where the CM energy is much larger than the rest mass. 

2.5 i 1 1 1 1 1 1 1 1 1 1 




q l i i i i i i i i i 

0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

a* 

FIG. 1: The CM energy E cra for the near-horizon collision. The solid, dashed, and dotted curves 
denote the collisions of a particle plunging from the prograde ISCO with (a) a marginally bound 
particle with I = II, (b) a particle plunging from the retrograde ISCO, and (c) a marginally bound 
particle with 1 = 0, respectively. For clarity, E cm is multiplied by y/l — a* in the vertical axis. 



V. PARTICLE COLLISION OF A PARTICLE ORBITING THE ISCO 

In this section, we deviate somewhat from the original idea of BSW [l|. We consider 
the situation where a particle orbiting the ISCO collides with another particle on the ISCO 
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instead of a near-horizon collision. In this case, we cannot take the near-horizon limit 
beforehand. Although we do not expect a compact expression for the general case, we can 
obtain a simple formula for the near-extremal limit. Using Eq. ( 12.1 Op with Eqs. ( I2.1ip -( l2.14p 
and (I4.5I) - (I4.7I) and y = M/risco; we can estimate the CM energy for the near-extremal 
Kerr black hole. If particle 1 is a particle orbiting the prograde ISCO and particle 2 takes 
a subcritical orbit, we obtain 

E cm r J V2e 2 - l 2 

2m ~ 2V6 3 V4 yi^I • ^ • > 

It is quite intriguing that the dependence (1 — a^) _1//6 on the spin parameter here is quite 
different from that for the near-horizon collision (1 — a^) _1//4 . The CM energy can be 
arbitrarily high in the near-extremal limit a* — > 1. However, this needs to be distinguished 
from the BSW effect for the near-horizon collision of plunging particles. 

As in the near-horizon case, we consider the "on-ISCO" collisions of a particle orbiting 
the prograde ISCO with (a) a marginally bound particle with / = l L , (b) a particle plunging 
from the retrograde ISCO, and (c) a marginally bound particle with / = 0. In the near- 
extremal limit a* — > 1, we obtain 

E cm ^ (3 
2^~ ^T^f 

where the numerical factor (5 = i/2e 2 — h/ (2 1//6 3 1 / 4 ) is calculated to be 
^2(2 + V / 2)/(2 1/6 3 1 / 4 ) ~ 1.77, 4 ■ 2 1 / 3 /3 ~ 1.68, and v / 2/(2 1/6 3 1 / 4 ) ~ 0.957 for cases 
(a), (b), and (c), respectively. 

For the general values of a, in < a* < 1, we can calculate the CM energy using 
Eq. (I2TT0D with Eqs. (l2TTll -( l2TT4l) and (Q]) -(1PI) . r = ri SCO and y = M/r lsco . The result 
is summarized in Fig. |2} where E cm / (2mo) is multiplied by y/l — a% for clarity. In this 
figure, the solid, dashed, and dotted curves denote cases (a), (b), and (c), respectively. The 
CM energy for two particles, either of which is a particle orbiting the prograde ISCO, is 
always below the upper bound (13. 8p for marginally bound particles. We should note that 
the dependence of (£ , cm /(2m ))(l — al) 1 ^ on a, in < a, < 1 is very weak and hence the 
empirical formula E cm /(2m ) ~ 1/ ^/l — dl is a very good approximation within a factor of 
2 for all values of a* in < a* < 1. Thus, this formula provides the typical value for the CM 
energy for the on-ISCO collision of a particle which orbits the prograde ISCO with a generic 
subcritical particle. It should be noted that the CM energy for the on-ISCO collision is 
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always smaller than that for the near-horizon collision of a particle plunging from the ISCO 
in the near-extremal limit because of the different dependence on a*. 

If we use Thome's bound a* = 0.998 for the spin parameter, E cm /(2m ) is calculated to 
be 4.11, 3.91, and 2.43 for cases (a), (b), and (c), respectively. This means that a highly or 
moderately relativistic collision can naturally occur on the ISCO around a rapidly rotating 
black hole in an astrophysical context. 

2 i 1 1 1 1 1 1 1 1 1 1 
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a* 

FIG. 2: The CM energy E cm for the on-ISCO collision. The solid, dashed, and dotted curves 
denote the collisions of a particle orbiting the prograde ISCO with (a) a marginally bound particle 
with I = II, (b) a particle plunging from the retrograde ISCO, and (c) a marginally bound particle 
with I = 0, respectively. For clarity, E cm is multiplied by y/l — a% in the vertical axis. 



VI. CONCLUSION AND DISCUSSION 

We have investigated particle collisions near the horizon and on the ISCO around a Kerr 
black hole. We have derived a general explicit formula for the CM energy near the horizon 
in terms of the energies and the angular momenta of colliding two particles on the equatorial 
plane. We have confirmed that the obtained formula includes known formulas as its special 

20 



cases. We have explicitly shown that although the effective potential around a maximally 
rotating Kerr black hole has a zero and maximal point on the horizon, it does not correspond 
to a real circular orbit. Then, we have studied the near- horizon collision of particles, either of 
which plunges from the ISCO. We have shown that the BSW effect occurs for such a collision 
in the near-maximal rotation limit and that the maximum value for the CM energy is the 
same within a factor of 2 as the upper bound for the marginally bound particles for which the 
angular momentum must be fine-tuned. We have also investigated the collision of a particle 
orbiting the ISCO with another generic particle on the ISCO and found that it is also the 
case that one can obtain an arbitrarily high CM energy in the near- maximal rotation limit, 
although this energy is smaller than the value for the near-horizon collision in this limit. 
The result implies that the BSW effect, which was originally proposed for the marginally 
bound particles with the fine-tuned angular momentum, is astrophysically relevant since 
the fine-tuning is naturally realized for ISCO particles in the standard accretion disks and 
extreme mass-ratio inspirals. Although the CM energy is bounded if the spin parameter of 
the black hole is bounded in an astrophysical context, the collision can still be highly or 
moderately relativistic near the horizon and on the ISCO around a rapidly rotating black 
hole. 

The present result naively suggests the following scenario. A highly or moderately rela- 
tivistic collision often occurs near the horizon of a rapidly rotating black hole in the context 
of the accretion disks and the extreme mass-ratio binaries. For the standard accretion disk, 
gamma rays with energy of several GeVs can be produced inside and around the inner edge 
of the disk, if the CM energies of protons and ions collisions are eventually converted to 
photons. These photons can have much higher energy than usual thermal photons. 

What is more intriguing is the high- velocity collision of compact objects around a super- 
massive or inter mediate- mass black hole. Here, the compact objects will collide near the 
horizon or on the ISCO with a high "relativistic gamma factor" E cm /(2mo)- The result will 
strongly depend on the kinds of the compact objects and the value of the relativistic gamma 
factor. For example, if two neutron stars collide with a sufficiently high gamma, it will result 
in the gravitational collapse to a black hole (e.g. 22J). If two white dwarfs collide with a 
sufficiently high gamma, they might be smashed, destroyed, and scattered away because of 
the CM energy much greater than the binding energy of the white dwarfs. Thus, the col- 
lision of the compact objects around a rapidly rotating supermassive or intermediate-mass 
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black hole provides a unique laboratory for the relativistic collision of black holes, neutron 
stars, and white dwarfs. The interaction of the compact object with the fluid or plasma 
near the horizon or on the ISCO might also be a striking phenomenon. The details of all 
these processes would not be so simple and should be investigated not only by analytical 
arguments but also by numerical simulations, including numerical relativity and general 
relativistic hydrodynamics. 

Finally, we speculate that peculiar signals originating from the highly or moderately rel- 
ativistic collision of particles, fluids and compact objects around a rapidly rotating black 
hole might be detected by the direct observation of black holes by means of electromagnetic 
and/or gravitational waves and/or neutrinos. For example, if two black holes collide with 
a sufficiently high gamma, a considerable fraction (as large as 14 ± 3% for head-on colli- 



sion 



and 35 ± 5% for zoom- whirl collision [24j, |25|) of the CM energy can be radiated 
away through gravitational radiation. Of course, because of the strong redshift, we cannot 
immediately expect that the energetic radiation can directly reach us. However, the emis- 
sion peculiar to such relativistic collisions will be redshifted and might still be observed in 
electromagnetic and/or gravitational waves and/or neutrinos. In this respect, the on-ISCO 
collision might be more advantageous to observation than the near-horizon collision. To 
investigate what signals would be observed from a highly or moderately relativistic collision, 
numerical simulations will be very powerful. 
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